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Abstract 

(N 

We define a "renormalized energy" as an explicit functional on arbitrary point con- 
figurations of constant average density in the plane and on the real line. The definition 
is inspired by ideas of |SS1[ ISS3J . Roughly speaking, it is obtained by subtracting two 
leading terms from the Coulomb potential on a growing number of charges. The func- 
tional is expected to be a good measure of disorder of a configuration of points. We give 
certain formulas for its expectation for general stationary random point processes. For the 
random matrix /3-sine processes on the real line (0 = 1,2,4), and Ginibre point process 
and zeros of Gaussian analytic functions process in the plane, we compute the expectation 
i - 1 explicitly. Moreover, we prove that for these processes the variance of the renormalized 

energy vanishes, which shows concentration near the expected value. We also prove that 
the (3 — 2 sine process minimizes the renormalized energy in the class of determinantal 
point processes with translation invariant correlation kernels. 
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1 Introduction 



The aim of this paper is to introduce and compute a function, called the "renormalized 
energy", for some specific random point processes that arise in random matrix models, and 
in this way to associate to each of these processes a unique number, which is expected to 
measure its "disorder". 

Our "renormalized energy", that we denote W, is defined over configurations of points 
lying either on the real line or on the plane, as the limit as TV" — > oo of 



2 sin V 



N 



logiV in dimension 1, 



nv({a,}) = -^ £ log 

i^j,a i: aj e[0,7V] 

WW({ai}) = 2nN 2 ^ E N {ai - a,) + log 2l ^tj\2 in dimension 2 > 

where En is an explicit Eisenstein series, and r\ is the Dedekind Eta function. 

This definition is inspired by that of the "renormalized energy" , denoted W, introduced 
by Sandier and the second author in [SSlJ in the case of points in the plane and in |SS3j in 
the case of points on the real line. The definitions for W and W coincide when the point 
configuration has some periodicity (this is where our new definition originates), and in that 
case they amount to computing a sum of pairwise interactions 

^G(ai - aj) 

where aj are the points and G is a suitable logarithmic kernel (the Green's function on the 
underlying torus); however they are not identical in general. We will give more details on the 
connection in Section [2] 

In [SS3] it is shown that in dimension 1, W is bounded below and its minimum is achieved 
at the perfect lattice Z. In dimension 2, the situation is more complex; it is also shown in |SSlj 
that the minimum of W is achieved, but it is only conjectured that this minimum value is 
achieved at the perfect triangular lattice or "Abrikosov lattice" according to the terminology 
of the physics of superconductors (which was the first main motivation for W in [SSl] where 
it was introduced). This conjecture is supported by the result that, among configurations 
of points which form a perfect lattice (of fixed volume), the renormalized energy is minimal 
if and only if the lattice is the perfect triangular lattice, i.e. with 60° angles (that result is 
shown in |SSlj based on the use of modular functions and results in number theory) . 

It is thus natural to think of W or W as a way to measure the disorder of a configuration 
of points. With this outlook, in dimension 1 the lattice Z is the most ordered configuration 
of points with prescribed density 1, while in dimension 2, it is expected to be the triangular 
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lattice (which is better ordered than any other lattice, say the square lattice). In addition, due 
to its logarithmic nature, W has some nice scaling and additive properties, which we believe 
make it a very good object. A further motivation for choosing to study W as opposed to any 
other total pairwise interaction function is that, as seen in [SS2] and [SS3], W arises very 
naturally from the statistical mechanics of Coulomb or log gases, which contain as particular 
cases the Ginibre and GOE/GUE/GSE ensembles of random matrices. In [SSl] , W was 
introduced and derived in the context of the minimization of the Ginzburg-Landau model 
of superconductivity. In |SS2| it was derived as a sort of limiting interaction energy for two 
dimensional Coulomb gases, and similarly in |SS3j with log gases. These works are based on 
analysis and energy estimates (upper and lower bounds). Both the questions we pursue here 
and the methods we use are quite different: they aim at obtaining explicit formulas for specific 
random matrix models. In particular, this is a way to compute some interesting statistics 
over random matrix eigenvalues, as initiated by Dyson-Mehta [DMj. We will comment more 
on this just below. 

Let us now briefly introduce the notion of a random point process. 

A (simple) random point process is a probability measure on the set of all locally finite 
collections of (mutually distinct) points in the space, cf. e.g. |DVJj . It can also be viewed as 
a random measure of the form = J2 P eA &p-> with the points p distinct and A discrete. 

Random point processes are essentially characterized by their "/c-point correlation func- 
tions" pk(xi, ■ ■ ■ , scjfc), which give the probability densities of finding k points at the locations 
x\, . . . , x\~. We will normalize our processes so that the average number of points per unit 
volume is always 1, which is equivalent to pi(x) = 1. 

Perhaps the most famous random point process is the Poisson process, characterized by the 
facts that the number of points in disjoint subsets are independent, and the number of points 
in any finite volume subset of the space follows a Poisson distribution with parameter equal 
to the volume of the set with respect to a reference measure. An important class of point 
processes is that of determinantal point processes, see |Sol] , |Ly| , [Jo] . |Kon05] . [HKPV] . 
|So2| . [B] and references therein. That class is characterized by the fact that the /c-point 



correlation functions are given by symmetric minors of a (correlation) kernel, cf. Section 3.1 
It is easy to see that lmijv_ 5 . 00 Wn = +oo for the translation invariant Poisson process, 
which means that it is "too chaotic" from the point of view of the renormalized energy. 
Let us list the (stationary) processes for which we show that W provides more meaningful 
information. 

We will be interested in one dimension in the /3-sine processes {(3 = 1, 2, 4) which arise as 
the local limit of the law of eigenvalues in random matrix ensembles with orthogonal, unitary, 
and symplectic symmetry groups (they are determinantal for (3 = 2 and Pfaffian otherwise). 
In two dimensions we will examine the "Ginibre" point process, which is also a determinantal 
process arising as the local limit of the law of eigenvalues of matrices from the complex Ginibre 
ensemble (i.e. square matrices with complex Gaussian iid entries), for further reference, see 
\Fo\ lAGZt [Me] ; as well as the random process of zeros of random Gaussian analytic functions, 
often denoted GAF, whose description can be found in [HKPV]. 

As our processes are always translation invariant, the 2-point correlation function can 
always be written in the form p2(x, y) = 1 — T-z(x — y) for the 2-point cluster function T2 (we 
will come back to this notation in Sections 3.1 and [5]). 



The main results we obtain are the following: 
• For general stationary processes we identify sufficient conditions on the process and its 
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2-point correlation function p2 for the existence of lirn/v-»oo EWjv, and give an explicit 
formula in terms of P2 which is (up to constants) 

lim EWjv = / log\x\T2(x) dx, 

7V->oo J R d 
with d = 1 or 2 according to the dimension. R 

• We apply this formula to the specific point processes mentioned above. 

• For the specific point processes above, we explicitly compute the limit of the variance 
of Wat and obtain that it is 0. This implies that for such processes, Wjv concentrates 
around its expected value, and converges in probability to limEWV as N — > oo. 

• We prove that in the class of determinantal point processes with translation invariant 
kernels in dimensions 1 and 2, limjv_ s>00 EWV is minimized by the processes whose cor- 
relation kernel is the Fourier transform of the characteristic function of the ball. A 
complete physical interpretation of this fact seems to be missing. In dimension 1 the 
optimization gives the (3 = 2 sine process, which can be seen as a heuristic explanation 
of the fact that this process is the universal local limit of (3 = 2 random matrix en- 
sembles. Indeed, such a local limit has to be translation invariant and determinantal, 
and it is natural to expect that it also minimizes an appropriate energy functional. In 
other dimensions, point processes with such specific kernels called "Fermi-sphere point 
processes" appear in [TSZJ as higher-dimensional analogues of the sine process. 

For our set of specific processes, we thus show that we can attribute to the process a unique 
number, which we compute explicitly. Whenever these numbers are distinct this implies 
that the processes are mutually singular (as measures on the space of point configurations). 
Moreover, we check that the processes that are expected to have the highest level of "order" 
or rigidity indeed have a lower value of lini7v->oo EWV. For example for the /3-sine processes 
we find 

limjv^oo E% = 2 - 7 - log2 for = 1 
lim^oo EW N = 1-7 for = 2 

lim^oo EW N = I - log 2 - 7 for (3 = 4, 

where 7 is the Euler constant. These three numbers form a decreasing sequence, as can be 
expected. 

In two dimensions, we obtain that the Ginibre process has a higher lim7 V _ 5 . 00 Wat, hence 
less rigidity, than that of zeros of Gaussian analytic functions, in agreement with the results 
and philosophy of [GNPSj. 

The values of limEWV for the (3 = 1,2,4 sine processes above equal twice the ther- 
modynamic "energy per particle" or "mean internal energy" for the log-gas with infinitely 
many particles, as obtained by Dyson in 1962 [Dy, I. IX, III. VI]. This is not surprizing as 
our definition of Wat essentially coincides with that of Dyson |Dy[ I. VI] once the points on 
the interval [0, N] are identified with points on a circle of same length. Furthermore, one 
can use this fact to infer finer asymptotic properties of WV as follows (we are very grateful 



1 Interestingly enough, the sufficient conditions involve the equality J Rd log\x\T2(x) dx = 1 that we check 
for the above mentioned processes. We expect it to hold for general /3-sine processes but we do not see an a 
priori reason for that. It is also not clear to us whether this condition has a physical meaning. 
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to Peter Forrester for the idea). If, instead of considering growing windows, one approxi- 
mates the /3-sine processes by the circular ensembles of particles on the unit circle with joint 
probability density Ili<j \ z i ~ z if '> one observes that in rescaled variables, Wn is simply 
— 2N~ l 1°S \ z i ~ z j \ + logiV, and its characteristic function can be immediately obtained 
from Selberg's formula for the partition function: 

7( ft — 2it \ rfl -X- P N \ 

E exp(it WW cul n = N u { ^ ra , N) , Z(P) - 1 + 2 ' 



z(P) ' ^' (r(i + §))*' 

see e.g. [Fo, Section 4.7.2] for the formula for Z(/3). Using Stirling's formula it is not hard to 
see that for any (3 > and tEl 

^limEexp (W 1 / 2 (n^ rcular - u(/3))) = exp (-^^y^) 



with 



and the expression of u(/3)/2 through the partition function Z(f3) coincides with that of the 
thermodynamic energy per particle. By Levy's continuity theorem this implies the central 
limit theorem 



N- 



hm Prob < — N , — < s > = — = / e x /z dx, s G 



It is natural to conjecture that the same central limit theorem holds for Wn on general 
/5-sine processes (constructed in |VVj ) as well, in particular providing asymptotic values 
limEVW = and limJV • Var(WN) = v((3) for any fj > 0. At the moment we do not 

know how to prove it, although in the (3 = 1,2,4 cases one may hope to obtain a proof via 
controlling the asymptotics of moments of Wn through explicit formulas for the correlation 
functions. 

In the follow-up paper [DMJ , prompted by the wish to analyze experimental data, Dyson 
and Mehta looked for a way of approximating the energy per particle, in the case /3 = 1, by 
averaging a statistic of the form ]P F(ai,aj) over windows of increasing length in the random 
matrix spectrum. The concentration requirement of the asymptotically vanishing variance 
led them to the conclusion that with their choice of F, the statistic had to be corrected with 
further interaction terms |DMl Eq. (109)]. Our renormalized energy Wn seems to be a better 
solution to Dyson-Mehta's problem, although one should note that Dyson-Mehta's statistic 
is purely local, while Wn involves long range interaction (between leftmost and rightmost 
particles in [0, N]). 

The paper is organized as follows: In Section [2] we give the precise definitions of W, the 
context about it and its connection to Coulomb energy, from |SS1| ISS3] . This leads us to the 
definition of Wn- In Section [3] we compute limits of expectations of Wn- First on the real 
line and for general processes, then for the f3 = 1,2,4 sine processes, then in the plane and 
for general processes, finally for the explicit Ginibre and zeros of GAF processes. In Section 
[4] we find minimizers of lim Wn among determinantal processes with translation-invariant 
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correlation kernels. In Section[5]we compute the limit variance of Wjy for our specific processes 
(and show it is 0). In Section [6] we gather some miscellaneous computations: the effect on 
limEVW of superposition and decimation of processes and the computation of expectation 
of limEVW f° r the f3 = 2 discrete sine processes. 
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related to Section |4j A. B. was partially supported by NSF grant DMS-1056390, S. S. was 
supported by a EURYI award. We also thank the MSRI for support to attend a workshop 
where this work was initiated. 



2 Definitions of W and Wn 

The aim of this section is to present a definition of W for a configuration of points either 
on the line or in the plane, which is directly deduced from that of [SSI] ISS3j . but depends 
only on the data of the points. We do not attempt here to fully and rigorously connect the 
next definition with that of |SSH ISS3j (since we believe it presents some serious technical 
difficulties) but the link will be readily apparent. 

We start by recalling the definitions from |SSlj . but rather in the form presented in [SS2J. 

2.1 Definition of W(j) in the plane 

In [SSI], W was associated to a discrete set of points (with asymptotic density) m in the 
plane, via a vector field j: if j is a vector field in M 2 satisfying 

(2.1) curlj = 2tt{v — m), div j = 0, 

where m is a positive number (corresponding to the average point density) and v has the 
form 

(2.2) v = 5 P for some discrete set Act 2 , 

peA 

then for any function x we define 

(2-3) W(j,x) ■= Km ( - / Xlil 2 + 7rlog?7 Vx(p) 

1-M Y JR2\U peA B(p, v ) ~ 

Definition 2.1. Let m be a nonnegative number. Let j be a vector field in M 2 . We say j 



belongs to the admissible class A m if (2.1), (2.2) hold and 



(2.4) ^ ^j- is bounded by a constant independent of R > 1. 
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Definition 2.2. The "renormalized energy" W(j) relative to the family of balls (centered at 
the origin) {-Br}r>o (and the number m) in M 2 is defined for j E A m by 

(2.5) W(j) : =limsup ^^ , 

where xr denotes positive cutoff functions satisfying, for some constant C independent of R, 

(2.6) |Vxh| < C, Supp( XR ) C B R , X r(x) = 1 if d{x, (B R f) > 1. 

Note that by scaling, the density of points can be changed to any, and the rule for that 
change of scales is 

(2.7) W(j) = m(w(j')-^logm), 

where v has density m and j' = -^j(-^=) (hence the set v 1 has density 1). 

This function was first introduced in |SS1| and derived as a limiting interaction energy 
for vortices of Ginzburg-Landau configurations. Independently of Ginzburg-Landau it can 
be viewed as a Coulombian interaction energy for an infinite number of points in the plane, 
computed via a renormalization. Many of its properties are stated in |SSlJ. In |SS2j it 
was directly connected to 2D log gases. We will give more details and properties of W in 
Section 12.31 



2.2 Definition of W(j) on the line 

In [SS3] a one-dimensional analogue to the above definition is introduced for the study of ID 
log gases, we now present it. The renormalized energy between points in ID is obtained by 
"embedding" the real line in the plane and computing the renormalized energy in the plane, 
as defined in [SSI]. More specifically, we introduce the following definitions: 

M denotes the set of real numbers, but also the real line of the plane M? i.e. points of the 
form (x, 0) G M 2 . For the sake of clarity, we denote points in R by the letter x and points in 
the plane by z = (x, y). For a function x on we define its natural extension x to a function 
on M 2 by x(x,y) = x( x )- Ir denotes the interval [—R/2,R/2] in R. 

5r denotes the measure of length on R seen as embedded in M 2 , that is 

/ <f5m = / f{x, 0) dx 

for any smooth compactly supported test function tp in M 2 . This measure can be multiplied 
by bounded functions on the real- line. 

Definition 2.3. For any function x on M., and any function h in M 2 such that 

(2.8) -Ah = 2-k(v - m5 m ) in R 2 
where m is a nonnegative number and v has the form 

(2.9) v = ^^<5p for some discrete set of points o/M, 
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we denote j = — V^/i := (ffe/i, —d\h) and 



(2.10) W(J, X) ■= lim - / xbf + 7rlog?7 Vx(p) 

^2 y R 2 VUj)SAjB(p!7?) ^ 

where x is the natural extension of x- 

Definition 2.4. Letm be a nonnegative number. We say j = — V^/i belongs to the admissible 



class A m if (2.8), (|2.9|) hold and 



(2.11) is bounded by a constant independent of R. 

We use the notation xr f° r positive cutoff functions over K satisfying, for some constant 
C independent of R, 

(2.12) IVxbI < C, Supp(xk) C I*, = 1 if N < ^ - 1. 
Definition 2.5. T/ie renormalized energy W is defined, for j £ A m , by 

(2.13) := limsup Xr) . 

Note that while W in 2D can be viewed as a renormalized way of computing ||/i||^i( K 2), 
in ID it amounts rather to a renormalized computation of ||/*|| #1/2 (jo- 
in one dimension, the formula for change of scales is 

(2.14) W(j) = m(W(j')-wlogm). 
where j corresponds to density m and j' to m = 1. 

2.3 Background 

We recall here some properties and background from |SS1| ISS3] . 

- Since in the neighborhood of p £ A we have curl j = 2n8 p — m(x), div j = 0, we have 
near p the decomposition j(x) = V -1 log \x — p\ + f(x) where / is locally bounded, and 



it easily follows that the limits (2.3), (2.10) exists. It also follows that j belongs to L p 
for any p < 2. 



loc 



Because the number of points is infinite, the interaction over large balls needs to be 
normalized by the volume as in a thermodynamic limit, and thus W does not feel 
compact perturbations of the configuration of points. Even though the interactions are 
long-range, this is not difficult to justify rigorously. 

The cut-off function xr cannot simply be replaced by the characteristic function of Br 
because for every p G A 

lim W(j, 1 Br ) = +00, lim W(j, 1 Br ) = -00. 
R<\p\ R>\p\ 
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- In dimension 1, there is a unique h satisfying (2.8) and for which W(— V^/i) is hnite. 
Thus W amounts to a function of v only. 

- In dimension 2, there is not a unique j for which W(j) < oo and it would be tempting 
to define W as a function of v only and not j by minimizing over j, however we do not 
know how to prove that the resulting object has good properties such as measurability. 

An important question is that of characterizing the minimum and minimizers of W. For 
the case of dimension 1, it is proven in |SS3] that the value of W does not depend on the 
choice of xr satisfying (2.12), that W is Borel- measurable over L^ oc (IR 2 , M 2 ) for p < 2, and 
that 

(2.15) minW = -vrlog(2vr) 

is achieved for v = ^ pgZ 5 P . 

For the case of dimension 2, it is proven in Theorem 1 of [SSI] that the value of W does not 
depend on {xb r }r as long as it satisfies (2.6), that W is Borel-measurable over L^ oc (IR 2 , R 2 ), 
for p < 2 and that min_4 m W is achieved and finite. Moreover, there exists a minimizing 
sequence consisting of periodic vector fields. The value of the minimum of W is not known, it 
is conjectured (see |SSlj ) that it is asymptotically equal to the value at the perfect triangular 
lattice. This conjecture is supported by the fact that the triangular lattice can be proved to 
be the minimizer among lattice configurations of fixed volume. 

In addition to arising from the study of Ginzburg-Landau |SSlj , W is naturally connected, 
as shown in ISS2, SS3|, to Vandermonde factors such as e~ 



(2.16) w n (xi, ...,x n ) = - ^log|xj - xj\ 

v£j 8=1 

where Xi, . . . , x n G M rf with d = 1 or 2, for any potential V (for example V quadratic). To 
explain this, let us introduce some more notation. We set 

(2.17) I{p)= f -\og\x - y\dpL{x) dpL{y) + ! V(x)d^(x). 

Jm. d xR d JR d 

It is well known in potential theory (see [ST] ) that, provided V(x) — log \x\ — > +oo as |x| — > oo, 
/ has a unique minimizer among probability measures, called the equilibrium measure - let 
us denote it fj,Q. It is characterized by the fact that there exists a constant c such that 

V V 

(2.18) + — = c in the support of /j,q, and + — > c everywhere 

where for any /x, is the potential generated by (a, defined by 

(2.19) u^{x) = -l \og\x -y\dp,{y). 



We may then set £ = + \ — c where c is the constant in (2.18). It satisfies 



(2.20) 



C = in Supp(/i ) 

C > in R d \ Supp(/io) 



It is easy to check that £ grows like V at infinity. The connection between w n and W is given 
by the following exact "splitting formula" , valid in both 1 and 2 dimensions. 
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Lemma 2.6 ( |SS3j ISS2] ). For any x\, . . . ,x n S d = 1 or d = 2, £/ie following holds 

n 1 n 

(2.21) w n (xi, . . . , x n ) = ™ 2 ^(W)) - " logn + -W(-V ± J ff / , 1 H „) + 2n V C(a 

r/ 7r ^ — * 



d TT 

i=l 



where W is defined in (2.3) or (2.10) respectively, and where 

(2.22) H = -2irA- 1 -n/zo J in M 2 , H'(n l/d x) = H( 



\i=l / 

where the equation is solved in M 2 and /iq is naturally extended into a measure on M 2 if d = 1. 
(Note that A -1 is the convolution with 1/ (27r) log | • \.) 

Sketch of the proof. We start by writing 

w n (x\, . . . ,x n ) = / — log \x — y\ dv{x) dv(y) + n I V(x)dv(x) 



where A denotes the diagonal of (M d ) 2 and u = X^ILi 'W The idea is to compute the right- 
hand side by splitting v as n/xo + v — n\i§. This way, using the fact that //o(A) = 0, we 
obtain 

w n (xi, . . . ,x n ) = n 2 IQu ) + 2n / U fl0 (x) d{y - nfi )(x) + n / V{x) d{y — nna){x) 



+ / -log|x - y\d{v - n /i ) {x)d{v - n /x ) (y). 
Since + \ = c + ( and since v and n//o have same mass n, we have 

2n J U>*(x)d(v -niM))(x) + n J V(x) d(u - nfi )(x) = 2n J (d(u-n ) = 2n J (dv, 
where we used that ( = on the support of hq. Thus 

(2.23) w(xi, . . . ,x n ) = n 2 I(no) + 2n J(dv + J - log \x - y\ d{v - n/j, Q )(x) d(y - nu )(y). 

We then claim that the last term in the right-hand side is equal to ^W(— V^-ff, ljgd). We 
define Hi(x) := H(x) + log \x — Xi\. We have Hi = — log*(z^ — n^o), with Vi = v — 5 Xi , and 
near Xi, Hi is C . It is then not difficult to deduce 

(2.24) J - log \x - y\d(u - n/j, )(x)d(u - nfx )(y) = Y^Hj(xi) - n J H(x)d/i (x). 
On the other hand, by definition ( |2.3[ ), ( |2.10 ) and using Green's formula, we have 

UiB(xi,r,) Z7T i JdB(xi,n) 0n 1 JM 2 \Ll i B(xi, V ) 

Using the decomposition H = Hi — log \ x — X{\ near each Xi, adding nlogry and then letting 
n — > 0, we arrive at the same result as the right-hand side of (2.24). This establishes the 
claim, and the final result then follows from the change of scales x' = n^l^x in W. 

□ 
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2.4 Calculation for points on a torus in dimension two 

The reason why we need another definition here is that, as we already mentioned, we wish 
to define W based on the knowledge of u, i.e. a set of points, alone. For a given v, there 



is no uniqueness of the j satisfying (2.1) or (2.8) (in fact the indetermination is exactly a 



constant, see [S.S2, Lemma 1.4]), and this makes it problematic to define W for points only 
in a measurable manner. However, if the configuration of points is periodic, then W can be 
computed nonambiguously from v only. 

The formula for this in dimension 2 is given in [SSlJ. If the periodicity is with respect to 
a torus T, and the number of points in that torus is n (denote them oi, . . . ,a n ) then there 
exists a unique (up to a constant) solution to 

(2-25) - AH {ai} = 2vr (j2 K - on T 

(i.e. a periodic solution) and we let 

(2-26) j{a i} = -V ± i? {tti} - 



It was proved in [SSlJ that W(j{ a .j) is the smallest of the W(J) over all j satisfying (2.1) 
which are T-periodic, and it was established (see formula (1.10)) that 



2 . 

where c is a constant, and G is the Green's function associated to the torus in the following 
way: 

(2.27) - AG = 2ir (d - ^) in T. 

In [SSlJ, (1.11), an explicit formula is given for this: 

i „2inp-(ai—aj) j 2i-Kp-x 

¥i P e(z2)*\{o} w x ■ • ■■■■ ! ' 

which can in turn be expressed using Eisenstein series. 

For simplicity we prefer to work with density 1 in a square torus, which is then necessarily 
TV := M 2 /(A^Z) 2 where n = N 2 . Compared to |SS1] . this will change the constants in the 
formulae. So for the sake of clarity we will include below a complete proof of the following: 

Lemma 2.7. Let ax, ... ,a n be n = N 2 points on TV- Let ji ai } be the Tjy-periodic vector 
field associated through ( |2.26[ ) to the configuration of points {a{\, extended by periodicity to 
the whole plane. Then W(j^ a .y) as defined in (2.5) is a function of ai, . . . ,a n only, which is 
equal to 

7T 



(2.29) W N (at, ■■-,a n ) = j^Yj ^ " a ^ + 7tR 
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where G is the unique solution of 
(2.30) - AG = 2tt ( 5 



1 

AT2 



in T 



N i 



u>zi/i J T cfcc = 0; and R is a constant, equal to lim^^o (G(x) + log |x|). 

Proof. First, arguing as in |SSlj Proposition 3.1, we have 



(2.31) 



W(j {at} ) = 



lim 



\j{a z }\ 2 + Trnlogr/ , 



7 AT | 2 yT^\Uf =1 B(ai,r7) 

which means we reduce the computation in the plane to a computation on the torus. The 



next step is a renormalized computation a la |BBH] . as skeched in the proof of Lemma 2.6 



Solving for (2.25) via the Green function (2.30), in view of the translation invariance of 
the equation, we can choose 

n 

(2.32) H {a . } (x) = J2G(x- ai ). 

i=i 

Let us denote 

(2.33) R(x) = G(x) + ]og\x\, 

it is well-known that R(x) is a continuous function. We denote R := R(0). 
Inserting (2.26) into (2.31), we have 

1 1 



(2.34) 



W (Hai}) = t^t 



lim 

r?— s-0 2 



T N \U? =1 B(ai,T,) 



W H {a t }\ 2 +vrnlog?? 



Using Green's formula for integration by parts, we have 



/ 



\VH. 



T N \Uf =l B(ai,v) 



{ai}\ 



[ (-&H {ai} )H {ai} + J2 I H {ai 



dH 



dn 



Inserting (2.25) and (2.32), we have 



/ 



\VH 



T JV \U^ 1 B(a i ,J7) 



{ad I 



T Ar \U™ =1 B(a i ,r / ) 



2ttG(x — aj) dx 



n n „ 



8H. 



{a;} 



<9n 



with n the inner unit normal to the circles. Using for the first term the fact that G is chosen 
to have average zero, and for the second term splitting G(x — aj) as — log \ x — aj | + R(x — aj), 



dH 



using the continuity of R(x) and lim^o f 9B ( a .^ —g^ 
we find 



2tt (from ( |2.25| )), and letting rj — > 



lim - 

r?— 5-0 2 



|VF( a j| 2 + irnlogn = n Y]i?(0) + ir Ya- log |aj - oy| + i?(aj - o 3 -)). 

T iv \U7 =1 B(a I ,r,) i=1 <?y 



Inserting into (2.34) yields (2.29). 



□ 
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It turns out that G can be expressed through an Eisenstein series. The Eisenstein series 
with parameters r G C and u,»£f (cf. [LaJ) is defined by 



(2.35) 



2in(mu+nv) 



Im[r) 



(m,n)eZ 2 \{0} 

Let us now define, for x G C, 

( 2 -36) ^at(x) := E Re ^ x / N ^ Im ^ x / N )(i) 



mr + n 



2 ' 



As in |SS1] we will also need another classical modular function: the Dedekind n function. 
It is defined, for r in the upper-half of the complex plane by 



(2.37) 



r,( r ) = q 1 / 24 f[ (1 - q k ) where q = e 



H-KT 



k = l 



The following holds 



Proposition 2.8. Let eti, . . . , a n ben = N points on Tjv, and let Wn be defined as in (2.29 ). 



En being defined in (2.36), we have 

W N {ai,...,a n ) = ^^^(ai-Ojj+ilog^ -2vrlogr?(z) 



2.38 



where points in the plane are identified with complex numbers. 



Proof. As in [SSlJ, G, introduced in Lemma 2.7, can be computed from (2.30) using Fourier 
series and this yields: 



(2.39) 



G(x)= J2 



e n 



p-x 



P& 2 \{0} 



2n\p\ 



We recognize here an Eisenstein series: 
(2.40) 



Note that the fact that G has zero average implies that 



(2.41) 



/ Ejf(x) dx = 0. 



There remains to compute the value of the constant R = lim x _-.o (G(x) + log |x|). For that 
we use the "second Kronecker limit formula" (see |La| ) that asserts that 

E u , v (r) = -27rlog\f (u-VT,T)q v2 /% 
where q = e 2l7TT , p = e 2l7TZ , z = u — vt, and 

(2.42) f(z,r) = gV12(pl/2 _ ^-1/2) JJ (1 _ q k p){1 _ q k /p) _ 

k>l 
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We use this formula with r = i, u = Re(^), v = Im(jj). In that case q = e , z 
and (where the bar denotes complex conjugation), and the formula yields 



X 

N 



(2.43) E N (x) 
As x — > we have p — > 1 and 



-27rlog 



pl/2_ p -l/2 



2k- 



iV 



and also (see ( 2.37|) ) 

q i/u n(! - 9*p)a - «vp) ~->o ? i/12 na - « fc ) 2 = ^ 



fc>i 



So as x — > 0, we have 

En(x) ~ — 27rlog 



rj(i) 2 2iir 



-2ir log 



27r| 
N 



4tt log rj(i). 



Combining with (2.40), it follows that R = —log ^ — 21ogr/(i). Inserting this and (2.40) into 
(12. 291) we get the result. □ 



We emphasize here that this formula is the exact value for W as defined in (2.5) provided 
we assume full periodicity with respect to Tjy. 

2.5 Calculation for points on a torus in dimension one 



We do here the analogue in dimension 1, i.e. we compute W given by (2.13) assuming that 
the point configuration is periodic with respect to T/v := M/(iVZ). We assume that there are 
n = N points a\, . . . , a at in TV, hence W is computed with respect to the average density 1, 
i.e. m = 2ir. There exists a unique (up to a constant) HfaA satisfying 

-AH {ai} = 2K (j2S {ai>0) -S^j 

in M? and which is iV-periodic with respect to the first variable, i.e H^ aj j(x + N,y) = 
H { ai }{x,y) for every (x,y) £ R 2 . We then set j{ a .y = -V ± H {ai y. Then, as in [SS3] we 
have 

Lemma 2.9. Let a\,...,a^ be N points on TV = M/(iVZ). Let j{ ai } be as above. Then 
W(js a .y) as defined in ( |2. 13 ) is a function o/ai, . . . , a/v only, which is equal to 



(2.44) 



W N (ai, . . .,a N ) = j^^^G(ai - aj,0) + ttR 



where G(z) is a restriction of the Green function o/Tjy x R, more precisely the solution to 



(2.45) 



A z G(z) = 2ir 5 ( o , 



so) 



N 



Z £ TV X 



with fj N G(x, 0) dx = 0, and R is a constant, equal to 

lim (G(x,0) + log\x\) . 

x— >0 

Here 5j N is the distribution defined over IV x R by J 5j N 4> = J T <f)(x, 0) dx. 
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Proof. The proof is analogous to that of Lemma 2.7 First we have 
(2.46) 



lim 



We again write 
(2.47) 



T N \ U-K) 2 y( Tw xK)\u£ 1 B((ai,0), J7 ) 



1=1 



VH M \ Z +TTNlog7 ] . 



with G(z) = — log \ z\ + R(z), R a continuous function. Using Green's formula we have 



\VH. 



(T N xM.)\ufL 1 B((a i ,0), V ) 



{ai}\ 



H {ai}8l N 



N 



+E 

i=i 



if 



dif. 



{li} 



dn 



From (2.47) we have 



lim 

7 ?- +Q ^(T iV xK)\U^ 1 B((a i ,0), 7? ) 



TV 

7=1 ^ T iV 



by choice of G. The rest of the proof is exactly as in Lemma 2.7, inserting the splitting of G 
and p!5l □ 



On the other hand we shall see below that we have the following explicit formula for G 
restricted to the real axis: 



(2.48) 



G{x,0) = - log 



2 sin 



7TX 



Note that a consequence of this and the zero average condition on G is that 



(2.49) 



/ log 


. TTV 


2 sin — 


lr N 


N 



dv = 0. 



The previous lemma and (2.48) lead us to the following 



Proposition 2.10. Let a±, . .. ,on be N points on TTV = M/(iVZ) and let Wn be defined as 

ir(ai - dj) 



in (2.44). We have 



(2.50) 



W N (ai,...,a N ) = --^^lo. 



2 sin ■ 



N 



IT log 



2vr 
iV' 



Proof. The first step is to prove ( 2.48 ) . This is done by solving ( 2.45 ) in Fourier series/transform. 
We choose the following normalization for Fourier transforms and series: 

7(0= / f(x)e- 2 ^<dx 
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Ck(f) = / f(x)e 2i n Lx dx. 



Then the Fourier inversion formula is f(x) = J f(£)e 2lnx '£ d£ and respectively 



2 ; 7r /.' 
N J 



Since G is HV-periodic in its first variable, we may define its Fourier transform as a Fourier 
series in the first variable and a Fourier transform in the second, i.e., for m £ Z and £ € M, 

G(m,0= [ [ G{x,y)e-^ L e- 2i ^dxdy. 



If G solves (2.45) then by Fourier transform, G has to satisfy 



^ 2 [^+e)(^" l .o = ^s i(W] -^ l ) J 



2tt- 

(0,0) - ]y 6 *.v 



It is direct to establish that 5j N = NS^ with 5^ by definition equal to 1 if m = and 
otherwise. Combining these facts, we obtain 



G(m,0 



/or (m,£)^(0,0). 



The undetermination of G at (0, 0) corresponds to the fact that G is only determined by 



(2.45) up to a constant. By Fourier inversion, it follows that 



1 ~ c 2i7rga:+2i7ri/g 



1 /• e 2i7r^a;+2i7r^ 

d£ + c = — > / 5 <i£ + c 



2vr(^ + a 



Using the formula / °° ^Sr^-dx = ne 2 z^ (cf. |PBMlj ) with b = 2ny and z = |m|/iV, we 
arrive at 

(2.51) 



cos(27rS3;) „ i I |m 
G(x,y) = V — ^e"^^ +c. 



m=l 



We next particularize to the case y = 0, and use Clausen's formula (cf. [Let Chap. 4]) 

cos(£;:e) 



(2.52) 
and thus we find 



E 

fc=i 



A; 



log 



2 sin 



for < x < 27T, 



G(x,0) = -log 



2 sin 



7TX 



+ C. 



The constant c should be chosen so that G(x, 0) dx = 0, which imposes c = 0, in view of 



the relation f T log |2sin ^| dx = 0, a direct consequence of (2.52). This establishes (2.48). 



In addition the value of R follows: R = linr^o G(x, 0) + log \ x\ = — log and inserting into 



(2.44), the result follows. 



□ 
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2.6 New definition for a general point configuration in the plane/line 

In order to define W for a general configuration of points, without referring to the current 
they generate, we choose to base ourselves on the formulas (2.38) and ( |2.50[ ) of the previous 
subsection. 

To get rid of some constants, we use however a different normalization^] and, for a given 
family of points {ai}, we define 



(2.53) W N = -~ £ 

i/i,Oi,aj6[0,AT] 



log 



2 sin 



71" (a. 



N 



+ logiV 



in dimension 1 



respectively 
(2.54) w N 



2irN 2 



y~] E N (ai - + log 



N 



i^j,a,i,aj£[0,N] 2 



2TTT](i) 



in dimension 2. 



The results of the previous sections suggest trying to define W as the limit as A~ — > oo of 
these quantities. More precisely, for a given random point process, Wn becomes a random 
variable, and we try to see whether it has a limit as N — >• oo. Again, we do not claim 
a complete rigorous connection between such a quantity (which only depends on the point 
configuration) and the original W, which we recall, was defined via a vector field j. We also 
emphasize that in (2.53) and (2.54) the number of points in [0, TV] resp. [0, A] 2 is no longer 
necessarily equal to A" resp. N 2 *. 

Let us comment a little on the minimization of Wn- Since our definition has been relaxed, 
the question of whether Wn achieves a minimum becomes unclear. However, in dimension 
1, observing that still VW(Z) = 0, we have the following statement (in dimension 2 no such 
result is available): 



Lemma 2.11. Assume ai,...,ap- are k points in [0, N]. Then 



(2.55) 



1 

AT 



#3 



2 sin 



ir(ai 



N 



+ log AT > 1 



k 
N 



log AT + -log-. 



Thus, for any point configuration Wn > (l — jf) log A" + ^ log jr, where k is the number of 
points in [0, N] . 

Proof. The proof is a simple adaptation of the proof of minimality of the perfect lattice in 
|SS3j . Let ai, . . . , cifc G [0, A/"], and assume ai < • • • < a^. Let us also denote u\ : i = cij+i — a^, 
with the convention afc + i = ai + N . We have X)i=i u i>i = Similarly, let = aj+ p — a«, 
with the convention a^+j = ai + N. We have Xh=i 



view the points ai as living on the circle M./(NZ). When adding the terms in ai — aj in the 



pN. By periodicity of sin, we may 
addii 

sum of (2.55), we can split it according to the difference of p = j — i but modulo N. This 
way, there remains 

1 _ . 7r(a,; — aj) . 2 vx - ^, L . ttm 



(2.56) 



A r 



E 1 ^ 



2 sin 



A r 



log AT = -— log 1 2 sin 



■p,i 



p=l i=l 



N 



+ log A", 



the choice of normalization in dimension 1 is based on (2.15 I 
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where [•] denotes the integer part. But the function log |2sinx| is stricly concave on [0, ir]. It 
follows that 



1 k l 



iru 



sin ■ 



i=i 



N 



< log 



2 sin 



log 



2 sin 



p-ir 



Inserting into (2.56) we obtain 

( 2 - 57 ) -^E 1 ^ 



2 sin 



7r(aj - a,j) 



N 



+ log N > 



2k 

"N 



[k/2] 

£ic, 

p=i 



2 sin 



P7T 



+ log N. 



On the other hand, we know that Wjv(Z) = which means that — j| ^l°g |2sin ^| + 



log iV = 0, but also, since this is true for arbitrary integers N, 

pir 



(2.58) 



[k/2] 

2 ^ log 

P =i 



2 sin 



k 



+ logk = 0. 



Inserting into (2.57) we are led to 
(2.59) 



2 sin 



ir(ai - dj) 



N 



k 



+ log N > -— log k + log N = 1 - — log N + — log — . 



N 



k 



N 



k . N 



N k 



□ 



3 Expectation of WW 

3.1 Expectation and 2-point correlation functions 



We now turn to evaluating Wat for random point processes. In view of its form (2.53)-(2.54), 
the expectation of the random variable WV can be computed from the sole knowledge of the 
second correlation function of the process. Indeed, recall that for any k > 1, the Appoint 
correlation function of a random point process in M. d is characterized by the property that 

(3.1) E V f ,„„..„, i)= / F(xi ,... Ml ,..^...^ 

— J(K d ) k 
i\,...,ik pairwise distinct \ > 

where the expectation is with respect to our measure on locally finite subsets X = {xj} C M. d , 
and F ranges over a suitable space of test functions, see e.g. j DVJ] . 

We note here that determinantal processes are a particular class of processes characterized 
by the fact that the correlation functions can be expressed as 

(3.2) p k (xi, ...,x n ) = det (K(xi,Xj)) iJe[lM 

for some kernel K(x,y), see |Solj . |Ly| , [Jo] , |Kon05j . |HKPVj . |So2j . [B] and references 
therein. This will be used later. 

Here we need specifically the two-point correlation function. In addition, for our function 
the formula simplifies when the process is assumed to be stationary (i.e. translation invariant). 
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From now on we make the basic assumption that we are dealing with a translation invariant 
multiplicity-free random point process in i or i 2 , of density 1 (i.e p\ = 1) with absolutely 
continuous correlation functions (hence, (3.1) holds). If p2(x,y) is the two-point correlation 
function of such a process, it is of the form r^ix — y) since the process is stationary. It is more 
convenient to work with the "second cluster function" T2 = 1 — r2 (we will give a general 
definition of the cluster functions in Section [5]) . 
Our basic assumptions thus imply 

(3.3) pi = 1, P2(%, y) = 1 — T%(x — y) for some function Ti- 

By definition of P2, the expectation of Wjv is, in dimension 1, 

vr(x - y) 



(3.4) 



1 

N 



log 



[0,iV] 2 



2 sin 



N 



T%{x — y) dx dy + log N 



(where we have used (2.53) and ( 2.49| )) and respectively in dimension 2 



(3.5) 



1 

'2ttN 2 



[0,N] 2 x[0,N] 2 



En{x - y)T 2 {x - y) dx dy + log 



N 



(where we have used (2.54) and (2.41)). The question is then whether these quantities have 
a limit as N — > 00. As we show below, this will only be true under additional assumptions 
which in particular ensure a sufficient decay of T2 . 

Finally, it would be most interesting to find natural conditions on the random points 
behavior (their spacing etc) that would guarantee the existence of a limit to EWat. 



3.2 Expectation for one-dimensional processes: theoretical formula 

Theorem 1. Consider a random point process X on the real line with the one-point correla- 
tion function pi(x) and a two-point correlation function p2(x,y), satisfying (3.3). Under the 
following assumptions 

1) sup„ eM |T a («)| < 00; 

2) there exist a sequence {o.n}n>i such that logN -C -C N l / 2 ~ e as N — )■ 00 (for some 
e > 0) and uniformly in A £ [a at, N — ajy] we have 



T 2 (v)log 



a N 



2 sin 



TTV 

iV 



dv = o(l) as N 



00; 



the following holds: 
-if 



if 



/OO l>OC 
T2{v) log \v I dv < 00 and / ^(u) dx = c 7^ 1, then EWtv — > 00 as N — >■ 00 
-00 J —00 

poo r&N 



T2M dx = 1 and 1 



T2(v) dv = o((logN) x ) for {aAr}Ar>i as above, then 



limTv^oo EWat exists and is finite if and only if J^ ^ T2(v) log \ v\dv converges, and if so then 



(3.6) 



lim EW N = log 2tt + 

JV-s>oo 



log |u|T2(f ) dv. 



We believe that Wjv should be bounded below or at least that the value —00 is in fact not taken. 
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Figure 1: Splitting of the domain of integration 



Remark 3.1. Condition 2) is satisfied by the stronger one: 

1 



(3.7) 



\T2(v )| dv = o 



log B 



as B —7- +oo. 



To see this it suffices to observe that on [atjsr, N — a at] we have | sin ^| > | sin ^jf-\ 



Proof of Theorem^ In view of (3.4) we need to compute 



(3.8) 



lim — 



[0,iV] 2 



log 



2 sin 



vr(x - y) 



N 



T2(x — y) dx dy + log N. 



For (x, y) G [0, N] 2 we denote u = x + y and v = x — y. We then split [0, N] 2 into the disjoint 
union of the following domains, see Figure 3.2, where cxn is in condition 2): 





= {(*> 


y) G [o,Af , 


K' ^ a A r ) Q A r ^ u ^ 2iV — Oat} 


Dx 


= {(x, 


y) e [o,Af , 


v > N — a at} 


Dy 


= {(x 


,y)e[o,iV] 2 


,v < -N + a N } 


D 2 


= {(x, 


v) e [o,m 2 , 


u > 2N -a N } 


D 2 , 


= {(x 


,y)e[o,iV] 2 


,u < a N } 


D 3 


= {(x, 


v) e [o,m 2 , 


Q-N < v < N — on} 


D 3 > 


= {(x 


,y) e [o,Af 


, — N + on < v < —on} 



We evaluate the integral in (3.8) over each of these domains successively. We start with 
the contribution of D[ . Making the change of variables a = Jr , b = ^jpt we have 



D' 



log 



2 sin 



7T X 



v) 



N 



T2 (x — y) dx dy 



a>0,b>0,a+b<°^f 



log |2 sin ir(a + b)\T 2 (N(a + b - 1)) da db. 
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Using assumption 1) and noting that in D\, sin7r(a + b) = 7r(a + b) + 0((a + 6) 3 ), and thus 
| log |2 siiivr(a + 6)|| = | log \2ir(a + b)\ \ + 0(^), we deduce 

7t(x - y) 



log 



2 sin ■ 



N 



T 2 (x — y) dx dy 
<iV 2 sup|T 2 | I 

J a 



log |27r(a + da db + O I ^ 

>0,b>0,a+b<^f V N 



a 



Using J Q r log sds = 0(s\ logs|), it follows that 



(3.9) 



/ log 


■ t(sg-v) 
2sm ^ 







^2(3; - y) cteciy 



a 7 



<ca^su P |t 2 |^| 



log 



N 



a; 



using a^v *C A^ 1 / 2 e . The estimate for the domains Z?i, D 2 , D 2 > is similar. For the contribution 
over D3, using the change of variables (x, y) — > (u, v) we have 



log- 



^a 



2 sin 



vr(x - y) 



T 2 (a;-y)dxdy = -(2iV-2a J v) 



N—ar 



T 2 (v)log 



2 sin ■ 



TTV 

iV 



This is o(N) when assumption 2) holds. The estimate on Dy is completely analogous. We 
have thus found that all contributions over D%, Di>, D2, Dy, D%, Dy are negligible. The be- 
havior of the integral will thus be determined by the contribution of Dq . 
Changing again the variables (x,y) into (u,v), we have 



log 



D 



2 sin 



ir(x - y) 



N 



1 



T 2 {x - y) dxdy = ~{2N - 2a N ) \ T 2 (v) log 



2 sin 



TTV 



But in D we have sin ^ = y(l + 0(§|)) as N -»• 00, hence log|sin^| = log | ^| + 0(^] 



iV 

2irv 



dv. 



Therefore, 



log 



D 



2 sin 



ir(x - y) 



N 



T 2 (x - y) dxdy 



, i-oln rax \ 

(N-a N )(-logN T 2 (v)dv + log \2irv\T 2 (v) dv J + sup |T 2 |0 



V 

A r 



Using oat € iV2, we easily deduce that if J"_^ ?2(f ) dv ^ 1 and log |27n;|T2(?;) dv < oo, 
then 

vr(a; - y) 



2 sin 



JV 



^2(2 — y] dx dy + log A" — > 00 



as A^ — )• 00, and we conclude as desired. If T 2 (y) dv = 1, then we may proceed and find 



(3.10) 



/ log 


o • Tr(x-y) 
2sm ^ 


'Do 





T 2 (x — y) dx dy + log A" 



log |27rv|T2(z;) dt> 



+ logN{l- T 2 {v)dv^j + o(^-logN^j +0 



JV 

AT 



and we conclude (returning to ( 3.4 ) and ( 3.8 ) , and using the assumptions) that ( 3.6 ) holds. □ 
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3.3 Specific computations on the line 

In this subsection, we use the result from the previous theorem to compute explicit asymp- 
totic values for EVW for some well-known point processes, namely the homogeneous Poisson 
process, and the /3-sine processes with /3 = 1, 2, 4. (Recall that VW was defined in such a way 
that for the lattice Z we have Wat = 0.) 

The homogeneous Poisson process satisfies p\{x) = 1 and p 2 {x,y) = 1, hence T 2 = 0. We 
immediately deduce from (3.5) that EWjv — > +00. Hence, the Poisson process can be viewed 
as having the 'value of W equal to +00. 

The /3-sine processes for /3 = 1,2,4 arise in random matrices as the local limits for random 
matrix ensembles with orthogonal, Hermitian, and symplectic symmetries, see |Me[ IAGZ1 IFoj 
and references therein. These are stationary processes whose correlations can be computed 
as follows: introduce the kernels 



(3.11) K^(x,y) 

(3.12) KW{x,y) 

(3.13) K^(x,y) 



sin7r(x — y) 
7t(x - y) 
/ sin7r(:r — y) 

tt(x - y) 
1 r< x -v) sint 

v Wo ~r 

I sin27r(x — y) 

2vr(x - y) 
I r2n(as-y) 

\ 2^ Jo ~ 



dt 



sgn(x - y) 



sint 



dt 



d sin27r(x 
dx 2tt{x — y) 
sin 2tt(x — y) 



d sin7r(x — y) \ 
dx tt(x — y) 
sin7r(x — y) 
7t(x - y) 



for (3 



for = 1 



for (3 



2tt(x -y) J 



where all the indeterminacies 0/0 at x = y are resolved by continuity. 

The (3 = 2 sine process is a determinantal process with kernel (x,y), thus from (3.2), its 
2-point correlation function is given by 

(3.14) p 2 (x 1 ,x 2 ) = det(K^{x u x j )) 

The correlation functions for the (3 = 1,4 sine processes have the form 

qdet (K^\xi,Xj 



(3.15) 



Pni^li • ■ ■ 1 -En) 



n 



1,2,... 



ije[l,n] 

where qdet denotes the quaternion determinant, see e.g. [Fo, Section 6.1.1] for a definition. 
Alternatively, the right-hand side of (3.15 ) can be expressed as the Pfaffian of a closely related 
matrix, cf. |Fol Proposition 6.1.5]. Random point processes with correlation functions of such 
form are often called Pfaffian, see [HI Section 10] and references therein. 

The three processes above satisfy p\ = 1, and their second cluster functions can be easily 
seen to be given by 

2 

for (3 = 2 



(3.16) 
(3.17) 
(3.18) 



T, 



( 2 ), 



(i) 



(4) 



v) 
v) 
v) 



sin7rv 

TTV 

sin ttv 

TTV 

sin 2ttv 
2irv 



1 d sin-7rf 
7r dv TTV 
1 d sin27Tf 
2tt dv 2ttv 



sint 



dt 



7T 



sgn(v) 







for (3 = 1 



2irv 



sint 



dt 



for (3 = 4. 
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Proposition 3.2. The /3-sine processes for (3 = 1, 2, 4 satisfy the assumptions of Theorem^ 
and 

lim E14V =1-7 for f3 = 2 

lim EVtV = 2 - 7 - log 2 for /3 = 1 

lim EWjv = - - log 2 - 7 for /3 = A, 
N— >oo 2 

where 7 is t/ie Euler constant. 

Before stating the proof, we recall some integrals of classical functions that we will need. 
We state them without a proof and refer to |PBM11 |PBM2j . 

(319) /„ — * = 2 

+0 ° sinw , , 7r 



(3.20) / logwcfo = -7 

7^ — 





f + °° fsmv\ 2 , vr 

r+°° /siniA 2 7r 
(3.22) y (— logvdv = --( 7 + log2-l). 

These formulas can be found in jPBMlj , they are respectively (2.5.3.12), (2.6.32.3), (2.5.6.9), 
and (2.6.32.7). Finally we need a few more integrals that are based on 

. [ x sint ... 7r f + °° sint 
«bz(x) = / at, sz(x) = £>z(a;) — — = — / at 

.In t 2 ./„ t 



and 



These are 



Ei(cc) = / — di. 
V —00 * 



CUi) I y +z2 dv = - S mh(z)Ez(-z) 



+ °° VSi(v) , 7T 



(3.24) / = K ' dv = -Ei(-z) 







v 2 + z 2 2 



f+°° Si(v) smv n Ei(z) - Ei(-z) 

cf. respectively (2.6.4.11), (2.6.2.10) and (2.6.4.16) in |PBM2j . 
Proof Proposition\3. S\ We start with the simplest case. 



Case /3 = 2: First it is easy to check that assumptions (3.3) as well as assumption 1) of 
Theorem jlj are verified. Also, we have T%(v) = 0(\v\~ 2 ) as v — > 00, and 1 — J^ N T^ = 
0((aAr) _1 ) = o((log N)^ 1 ), thus assumption 2) is implied by (3.7). According to Theorem[l] 
we then have 



00 / ■ \ 2 

' sin TTV 



lim EW N = / log \2ttv\ dv. 

AT-j-oo \ TTV 
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From (3.22), we obtain 
(3.26) 



' sin TTV 



TTV 



log \2ttv\ dv = 1 — 7, 



hence the result. 

Case (3 = 1: Similarly to the case of /3 = 2, one checks that T2{v) = 0{\v\~ 2 ). Indeed, note 
that /;H -fei dt - f = j£f ate* dt = ^1 _ ;+« «t = o^,-!). Assumptions Q and 

1), 2) are then verified as in the (3 = 2 case. Next we check that J = 1. By evenness and 
integration by parts 



T^'(v)dv = 2 



sin 7ri> 



TTV 



1 5 



Sin 7Tf 



TT 5f 7Tf 



sin t , TT 



t 2 



oo / ■ \ 2 

' sin TTV 



TTV 



dv-l = l. 



The convergence is also fast enough, since I^u) = 0(\v\ 2 ). 
According to Theorem [T] we thus have 



lim EW N 

7V->oo 



log \ 2ttv\ 



oo 
oo 



log |27TU| 



sin TTV 



TTV 



sm7rw 



TTV 



1 d sin7ru 

TT dv TTV 



sin t tt 
t a* - -^sgn(v) ) \ dv 



Sill 7TD 



ld_ 

TT 8v \ TTV 



sin i , tt 

dt 

t 2 



dv. 



Integrating by parts, we are led to 
(3.27) 



lim E14V = 4 / log(2vrv) 

N—too 



sm TTV 



TTV 



dv + 



TT 



sm TTV 



TTV 



- 1 

2 

TT Jo 



nv sint tt\ du 
2 J V 



+oo 



SIIITTV 

log(27rw) dv. 



In view of (3.26) the first term on the right-hand side is equal to 2 — 27. The second is equal 
to 



+00 



sin TTV \ dv 
1 si[TTV) — = - 

TTV J V TT Jq 



sin u \ du 
1 si(u)- 



u 



With (3.23) and (3.24), 



2 f + °° 
n Jo 



— lim 

TT z->0 



sinvru \ dv 

1 Sl{TTV) — = lim hil{ — Z) 

TTV I V z->0 



si(u) sin it si(u)u 



u 2 + z z 



sinh(z) 



v? + z 2 



du. 



1 • 



On the other hand, near x = one has 



(3.28) 



Ei{x) = 7 + log |x| + 



n=l 



nn! 
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hence the above right-hand side limit is 0. The second term in (3.27) thus vanishes. From 
(3.20) and (3.19), the third one is equal to 7 — log2. This concludes the case (3 = 1. 
Case (3 = 4; Assumptions (3.3) and 1) are easily verified, we proceed to 2). We may write 



T2(v) 



'sin27riA 1 d sin27ru 
2ttv J 2-7T dv 2irv 



2lTV 



sini 



dt 



7T 



sgn{v) 



1 d sin27Tf 
4 dv 2ttv 



sgn(v). 



The first part is 0(\v\ 2 ) just as in the case above, so it remains to check that 



d sin27ry 
dv 2irv 



log 



a N 



2 sin 



TTV 



dv = o(l) 



uniformly in A £ [ajy, N — ajv]. Integrating by parts, we have 
' 1 d sin27rti 



QJV 



dv 2ttv 



log 



2 sin 



TTV 

iV 



sin 2-7T A 
2irA 



log 



2 sin 



dv 

7rA 



N 



sin 2naN 



27TU 



N 



log 



2 sin 



A sin 2irv cos ^ 

— dv 

aN 2v N sin ^ ' 



In [a at , A] we may bound from below sin ^ by sin ^jfi- , which is asymptotically equivalent to 



^jf- as N — > 00. Hence the integral on the right hand side may by bounded by ^- J a * < 



= o(l) in view of assumption 4). The other terms are also easily found to be o(l) by 
a similar argument. 

We next check that J T% = 1 with fast enough convergence. By evenness and integration 
by parts 



T^\v)dv = 2 




1 d sin2-7Tf 
2ir dv 2ttv 

1 sin27raAr 
2ir ajv 




27rcw sint 



dt = 1 + OUaNr 1 ). 



According to Theorem [T] we thus have 



lim EW N 



log \ 2ttv\ 



sin 2ttv 



1 <9 sin27ru 



2ttv J 2tt dv 2ttv 
Using evenness and integration by parts as above, we find 



2lTV 



sint 



dt dv. 



lim EW N 



sm27rtA 1 

log \2irv\ dv H — 

2irv ) 



7T 



sin 2ttv 
2ttv 2 



2ttv 



sint 



dt dv. 



By (3.22) the first term on the right-hand side is equal to —7 — log 2 + 1. By change of 

,Ei(z) - Ei(-z) 



variables, the second term is equal to 



Si(u) sin-u 



du 



— lim 

7T z->0 



+00 



Si(u) sinn 
u 2 + z 1 



du 



- lim e 

4 z-s>0 



by (3.25). Combining with (3.28) we find that the second term is equal to \ and we conclude 
the proof. □ 

As expected lmijv-s.oo IEWtv decreases as ft = 1,2, 4 increases. 
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3.4 Expectation for two-dimensional processes : theoretical formula 

In the plane, the computations are easier because we can take advantage of the fast (expo- 
nential) decay of the correlation kernels. 

Theorem 2. Consider a random point process X in the plane, with a one point correlation 



function pi(x) and a two-point correlation function p2(x,y), satisfying (3.3). We identify the 
plane with the complex plane C. Under the assumption 

/ \v\ k \T 2 {v)\dv < oo for k = 1,2, 3; 
Jr 2 

the following holds: 

- if I T 2 (v) dv = c 7^ 1, we have EWjv — >• 00 as N — )• 00; 
Jr 2 

-if Tzlv) dv = 1 and 1 — J,_ N N]2 T 2 {v) dv = o((log iV) -1 ), then limAr_ s . 00 EWV exists and 

Jr 2 ' 
is finite if and only if J R2 T 2 (v) log \ v\dv converges, and if so, then 



(3.29) 



lim EW N = / log\v\T 2 (v)dv. 

N^oo J R 2 



Proof. Returning to (3.5) we have to compute 



1 



lim 

AT->oo 2ttN 2 



[0,Af] 2 x[0,Af] 2 



E{x — y)T 2 (x — y) dx dy + log 



N 



2nr](i) 2 ' 



Making the change of variables (u, v) = (x + y, x — y), we have 
/ E N {x-y)T 2 {x-y)dxdy = \ [ [ 

J\0,N] 2 x\0,N} 2 4 Jve\-N,N] 2 Ju 



vG[-N,N] 2 Ju&S N (v) 



E(v)T 2 (v)dudv, 



where Sn(v) = {x + y : x £ [0, N], y G [0, N] 2 , x — y = v}. We may compute that |5jv(v) 
4iV 2 - 2N\vi\ - 2N\v 2 \ + \vi\\v 2 \, so 



(3.30) / E N {x - y)T 2 {x - y)dxdy 

l{0,N] 2 x[0,N} 2 

' 1 {AN 2 -2N\vi\-2N\v 2 \ + \vi\\v 2 \) E N (v)T 2 (v)dv. 

-N,N] 2 



Next we return to (2.43) where / is given by (2.42) and perform an asymptotic analysis 
as N — > 00. We have 



pl/2 _ p-V2 = e wf _ e -^f = 2iiT^ + O 



X 

N 



\x\ 

w 



while, since p = 1 + O(^) we may write (with q = e 



-2tt\ 



a 1 ' 12 nu - 9 fc p)(i - q k /p) = ( q 1/12 " J (1 + o(p - 1)) = v(i) 2 + o , 

k>l \ k>l J ^ ' 
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hence 



V !) 1 + \:v 



Inserting into (2.43) and combining with e n ^ Im N^ 2 = 1 + O(Jft) we obtain 



(3.31) 



En(x) = — 2-7T log \x\ — 2tt log 



2irr](i) 



N + °\N 



as N — > oo. 



Inserting this into (3.30), we are led to 



[0,N] 2 x[0,N] 2 



En(x - y)T 2 (x - y) dxdy 



-N,N] 2 

Therefore 
1 



1 



1 



1 



N 2 - -JVM - -N\v 2 \ + 7 |ui||u 2 | -27rlog|w| - 2vrlog 



2irr/(i 



N + °\N 



T 2 (v) dv. 



2vrJV 2 



[0,7V] 2 x[0,A r ] 2 

1 N 
log 



En(x - y)T 2 (x - y) dx dy + log 



N 



2-KT]{i) 2 
+ 



-N,N] 

Using the assumption we have 



■N,N] 



log \v\ + O 

2 \ \ 1\ 



,v 1 + lf ) + ° {W^ J ^ rU!r 



"I ^ +0 ^ ) ; / • ■ 



4 


-N,N] 2 


9Uv)dv 


= o(l), 


4 


-N,N] 2 


^T 2 {v)dv 


= o(l), 


4 


-N,N] 2 


^T 2 {v)dv 


= o(l), 



I[-N,N] 2 9 lQ g \ V \ T ^( V ) dv = ol 1 )' 



Jj_jvjv]2 % log |t>|T 2 («) dv = o(l). 



It follows that, as N — > oo, 
1 



2vriV 2 



[0,A r ] 2 x[0,A r ] 2 



log 



-Eiv(^ - y)T 2 (x - y) dxdy + log 



N 



2-KT)(if 



N 



2vrr ? (i) 2 

The result then easily follows. 



T 2 (v) dv) + log |u|T 2 (t;) dv + o(l). 

-7V,Af] 2 / ./[-Ar.jV] 2 



□ 
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3.5 Specific computations in the plane 

We turn to computing that limit for two specific processes. The first one is the determinantal 
random point process with correlation kernel 

(3.32) K(x,y) = e-^ 2 +\y\ 2 -^v). 

This process arises in random matrices as the local limit of the complex Ginibre ensemble, 
see e.g. |Fol Proposition 15.2.3], and thus it is sometimes called the Ginibre point process. 
From the determinantal structure of the correlation functions, cf. (3.2), we have p2(x,y) = 
1 — \K(x,y)\ 2 = 1 — e~' K \ x ~y\ and T%(v) = e _7r ' v ' . This easily satisfies all the assumptions of 
Theorem pi (in particular J T2 = 1) and we obtain 



Proposition 3.3. The determinantal process with kernel (3.32) satisfies 

lim EWn = -~ (7+ log7r) . 

This statement can be compared to a computation done in |Jaj . see also |Fo[ Ex.l5.3.1(iv)]. 

Proof. According to Theorem [2] it suffices to compute 

f _ I ,2 f°° _ 2 f°° 1 

/ log I v I e w ' v ' dv = logre * T 2-nr dr = I -(logs — log 7r)e s ds, 
Jr 2 Jo Jo 2 

using the change of variables s = irr 2 . We have J °° e~ s log sds = —7, and the result follows. 

□ 

The second one is the process of the zeros of a Gaussian analytic function (often denoted 
GAF). It consists of the random zeros of the analytic function Y^n=o wnen the £ n are 

i.i.d Gaussians suitably normalized, and it is a stationary process in the plane. The general 
background can be found e.g. in [HKPVJ. 

The second cluster function for the process, when the density pi is taken to be 1, is given 
according to [FH] by 



T 2 {x) = l-h(^f 
where 

1 d 2 

h{x) = 1 + - — (x 2 (cothx - 1)) . 

2 dx z 

It is easy to check that the assumptions of Theorem [2] are satisfied, and we deduce 
Proposition 3.4. The "zeros of Gaussian analytic functions" process satisfies 

lim E14V = - X (l + logvr). 

N^roo 2 
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Proof. To check that we may apply Theorem [2] we first compute 

|2- 



T2(v) dv 



1-h 



1 - h 



7T\V\ 



dv 



2nr dr 



(1 - h(u)) du 



d 

dx 



(x 2 (cothx — 1)) 



[2cc(cothx - 1) + x 2 (l - coth 2 x)]| 



where we have used the change of variables u = nr 2 /2 and the asymptotic relation cothx ~ - 
as x —7- 0. It is also easy to check that the convergence of J T2 is exponential, hence fast 
enough, and we may apply Theorem [2} 
This yields 



lim EWV 

N-*oo 



I log |u|T2(v) dv 
Jr 2 



2tt 
2 



log r I 1 — h 

f2u 



( irr 2 



r dr 



log 







V 2 

(1 - h(u))du 



2u\ d 2 o 
log — -r^\ u (cothn — lJJ du. 
it ) du 1 



Let us now compute 
log 



^ (u 2 (cothu — 1)) du 



7T 



du 2 



log — 1 — (x 2 (coth x — 1)) 
n J dx 



- 1 d , 2 



it dti 



(it (cothu — 1)) du 



2e Z" 00 
log— (2e(cothe - 1) + e 2 (l - coth 2 e)) - [x(cothx - l)]f - / (cothu-l)du 



log — (1 + O(e)) + 1 + O(e) - [logsinhx - x 

TT 



00 
e 



Taking the limit e — > 0, we conclude 



lim EHV = --(1 + logTr). 



□ 



It is well known (cf. [NSJ and references therein) that the "GAF process" is more "rigid" 
than the "Ginibre point process" (cf. also recent work [GNPSJ). We just demonstrated here 
it has more order via the renormalized energy. 
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4 Optimization over determinant al processes 



As explained in the introduction and Section 2.3 the question of minimizing W is an impor 



tant one, and open in dimension 2. It thus seems interesting to try to minimize limjv_ 5 . 00 EWjv 



as expressed in (3.6) and (3.29), over a subclass of processes. 

In this section we show that we can characterize the minimizer of this expression over the 
class of determinantal random point processes whose correlation kernel K(x, y) is Hermitian 
and translation invariant, i.e. K(x, y) = k(x — y) for some function k. For those processes, we 
have T2 = k 2 . Note, however, that the important determinantal process with kernel (3.32) is 



not in this class: while all its correlation functions are translation invariant, the correlation 
kernel is not. 

We prove the following statement. The proof relies on a rearrangement inequality. 

Theorem 3. Let 1C be the class of determinantal processes on the real line, respectively the 
plane, with self-adjoint translation-invariant kernels K(x,y) = k(x — y), and k(v) G L 2 (R d ), 
d = 1 or d = 2, such that 

1) p\{x) = k(0) = 1 and J Rd k 2 (x) dx = 1; 

2) J Rd log |x|/c 2 (x) dx < 00. 

LetT(k) = J^ d log \x\ k 2 (x) dx. Then for any process from K, with correlation kernel K (x , y) = 
k(x — y), we have 

where B is the ball centered at of volume one, and ^ is the Fourier transform. Thus, on 
the real line J-(k) is minimized over K, by the (3 = 2 sine process, while on the plane the 
minimizing determinantal process has the kernel given by k(v) = 1b(v) = , where 

J\ is a Bessel function of the first kind. 

Remark 4.1. Condition 1) says that our processes have density 1, which we have assumed 
throughout, and that limAr^ooEWV is finite, cf. Theorems^ and\^ The functional J- coin- 
cides with lini7v->oo EWiv ; given that the decay assumptions are satisfied. 

Remark 4.2. Numerical integration shows that in dimension 2, J- (1b) ~ —0.65, which is 
greater than linijv^oo EWat for the Ginibre ensemble (= — ^(7 + log7r) « —0.86) and for the 
zeroes of the Gaussian analytic function (= — ^(1 + log7r) —1.07). Thus, the two latter 
processes are "more rigid". 

Proof of Theorem^ Let us denote by / the inverse Fourier transform of k, and by Tr the 
integral operator corresponding to K via Tx(^p) = J K(x,y)(p(y) dy. We have 

t k(v) = k(x- y)ip(y) dy = k*p. 



By Macchi-Soshnikov's theorem, see [Solj . any self-adjoint translation-invariant correlation 
kernel of a determinantal process gives rise to an integral operator with spectrum between 
and 1. Hence, the spectrum of Tk is in [0,1]. Since Tk is also the convolution by k, 
this implies that f = k takes values in [0,1]. Moreover / G L 2 with J \f\ 2 = J k 2 = 1, 
and J f = k(0) = 1 from assumption 1). A function / with values in [0,1] which satisfies 
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f f = f f 2 = 1 can only be a characteristic function of a set A, denoted / = 1a, with A 
measurable of measure 1. Writing Ua for the corresponding k, we may write 

k A (x) = J e- 2i ^ x l A (0^. 

There remains to optimize over A, measurable set of measure 1, 

(4.1) 1(A) := / k A (x) log |x| da. 

Let A be measurable of measure 1, such that 1(A) < +oo. Noting that kA is an L°° (and also 
continuous) function, for every a > the integrals 



/ 



.2, ,1-Ix| 



A;^ (x) fix 

a 

then also converge, by comparison. Given any r£l, (e Th — l)/h converges to r monotonically 
as /i — >• (it suffices to check that this function is increasing in h). It then follows that in 

1 \x\~ a 

each of the domains \x\ < 1 and \x\ > 1, we have — '— ■ > log \ x\ as a — > 0, monotonically. 

Splitting the integrals as sums over these two regions, it follows by monotone convergence 
theorem that 

(4.2) 1(A) = lim / 4(x) 1-|x| - dx. 

We then remark (in view of the formula f * g = fg) that k\(x) is the Fourier transform of 



Ia(x) 



[ i A (y)iA(y-x)dy = \(A + x)nA\. 



We next claim that Ja is a continuous function. First, consider the case where A is an open 
set. Then, as x — > xq we have lyi(y)l^(y — »)—)■ l J 4(y)lyi(y — xq) almost everywhere, by 
openness of A, while — x)| < 1a(2/) and 1a € The claim is thus true by 

dominated convergence theorem. Second, if A is a general measurable set, by outer regularity 
of the measure we may approximate it by an open set U such that A C U and |E/\A| < e. 
Then it is immediate that for any x, \/a(x) — fu( x )\ = \\(A+x)nA\ — \(U+x)nU\\ < 2e. Since 
fu is continuous and e is arbitrary, it follows that /a has to be continuous too as a uniform 
limit of continuous functions. The claim is proved. We also note that /a(0) = |A| = 1. 
The next ingredient is that in dimension 1 

(4.3) (1 = <$ " 2T(-a + 1) sin ^(2vr|C|) Q - 1 , 

while in dimension 2 

(4-4) (1 HxT"«) = 5 - vr^^^l^r- 2 

For a reference see e.g. |Ed|. Chapter 5], or jSchl page 113]. 

Let us continue with the one-dimensional case. We deduce from the above facts that 

(4.5) I a (A):= [ kKx) 1 '^ dx=- [ / A (e)(«5-2r(-a + l)sin^(2vr|e|r- 1 ) d£. 

Jr a a J R \ 2 J 
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The relation can be justified by convoluting S — 2Y(—a + 1) sin ^p(27r|£|) a 1 with a Gaussian 
kernel approximating 5 at scale e, using the fact that J fg = J fg in the Schwartz class, the 
continuity of /a and then letting e — > on both sides. Moreover, this argument shows that 
/ /A(£)l£| a_1 d£ is convergent. 



Using /a(0) = 1 and Fubini's theorem, we may rewrite (4.5) as 

(4.6) I a {A) = - (l-2r(-a + l)sm^ / l A {y)l A {v ~ x)(2vr|x|) a - 1 dxdy) 

= I _ lr(-a + 1) sin ^ f l A (y)l A (z)(2Tr\y - z\) a ~ l dz dy. 
a a 2 J R 2 

Notice that 4r(— a+1) sin ~ | as a — )• so for a small enough, — fr(— a+1) sin ^(2^1^ — 
2|) a_1 is increasing in \y — z\. Now Riesz's rearrangement inequality (see |LL| Theorem 3.7]) 



asserts that a quantity of the form of the right-hand side of (4.6) is always decreased by 
changing 1 A into its symmetric rearrangement (1a)* = 3-A*- This means that for all a small 
enough, I a (A) > I a (A*). But 1(A) = lim^^o Ia(A) hence the same is true also for / i.e. 
1(A) > I (A*). The symmetric rearrangement A* of A is the ball centered at and of volume 
\A\ = 1. We have thus found that f k\(x) log |x| dx is minimal when A is the ball centered 
at and of volume 1. In dimension one, the Fourier transform of lri ii is mnnx which 

corresponds to the determinantal process with K(x,y) = b ™^^ , that is the sine process 
(for P = 2). 



In dimension 2, the argument is exactly parallel, starting again from (4.4). □ 



5 Computations of variance of Wn 

In Section [3] we dealt with the expectation of W. In this section we turn to examining its 
variance in the sense of computing limAr^oo Var(WV) for the same specific random point 
processes. 

In what follows we will need the formalism of (higher) cluster functions to efficiently deal 
with the /c-point correlation functions for k = 2,3, 4; we refer to |TWl IFa| for details and 
further references on this formalism. 

For any nonempty subset S = . . . , i^} of {1, . . . , N} we write ps = pk(xi 1 , . . . , Xi k ), 
where pk is the Appoint correlation function, and define the n-point cluster function as 

(5.1) T n (xi, ...,x n ) = ^(-l) n ~ m (m - l)\p Sl ... PSm 

with the sum running over all partitions of {1, . . . , n} into nonempty subsets Si, ... , S m . From 
the T n , the p n can be recovered through the reciprocal formula 

(5-2) p„ = ^(-ir- m T Sl ...T 5m . 

If a random point process is determinantal (cf. (3.2)) with correlation kernel K, then (see 
e.g. [TWJ) for any k > 1 

(5-3) T k (xi,...,x k ) = - ^2 K(x a{1) ,x ai2) ) ...K(x a{k) ,x a{1) ), 
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where S k denotes the symmetric group on k symbols. 

If the correlation functions of a point process are given by quaternion determinants (cf. 



(3.15)) then (see e.g. fFW\ ) for any k > 1 



(5-4) T k (x 1; ...,x k ) = ^-Tr ^ K(x a{1)l x ct(2) ) • • . K(x a{k) , x a{1) ). 



2k 

o-es fe 



Lemma 5.1. We have 



Var j ^ Gn^cli - aj) J = h -\ h h 



where 








(5.5) 


h = 


2 / Gn(x — y) 2 dx dy 

J[0,N] 2 




(5.6) 


h = 


-4 / G N (x - y)G N {x 

J[0,N] 3 


— z)T 2 {y, z) dx dy dz 


(5.7) 


h = 


2/ G N (x-y)G N (z- 

J[0,N] 4 


t)T2(x, z)T 2 (y, t) dx dy dz dt 


(5.8) 


h = 


-2 / G N (x - y) 2 T 2 (x 

J[0,N] 2 


y) dx dy 


(5.9) 


h = 


4 / G N (x-y)G N (x- 

J[0,N] 3 


z)Tz{x, y, z) dx dy dz 


(5.10) 


h = 


- / G N (x - y)G N (z - 

J[0,N] 4 


- t)T^(x, y, z, t) dx dy dz dt, 


where Gn(x) 


= - log 


2sin^ in dimension 1, 


resp. Gn(x) = t^En(x) defined in 



dimension 2. 

Proof. Expanding the square, we have 

( 5 - n ) ^ G N {di - aj) = ^2 G N(ai - aj)G N {a k - a,) 

\iy£j,a it a,j£[0,N] J i,j,k,l p.d. 

(5.12) + ^ G N (cii - aj)G N (ai - ai) + ^ G N (ai - aj)G N (a k - at) 

i,j,l p.d. i,j,k p.d. 

(5.13) + ^2 G N (ai - aj)G N (aj - ai) + G N (ai - aj)G N (a k - aj) 

i,j,l p.d. i,j,k p.d. 

(5.14) + ^2 G N{ai - aj) 2 + ^2 G N(ai - a j )G N (a j - a*), 

where the sums are still taken over points in [0, N], and p.d. stands for "pairwise distinct". 



Since Gn is even, it is clear that all the sums in (5.12) and (5.13) are equal, and the sums in 
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(5.14) as well. Using /c-point correlation functions (cf. (3.1)), we thus may write 

(5.15) E ( G ^(^ ~ a j) 

\i^j,ai,aj£[0,N] 



[0,7V] 4 



Gn(x — u)Gn(z — t)p4(x, y, z, t) dx dy dz dt 



+ 4/ G N (x -y)p 3 (x,y,z)dxdydz + 2 G N {x - y) 2 p 2 (x - y) dx dy. 

./[0,7V] 3 ./[0,7V] 2 



It is now convenient to express this in terms of the cluster functions T^, using (5.2), which 
yields 

p 2 (x,y) = T 1 (x)Tx(y) -T 2 (x,y) = 1 -T 2 (x,y) 
P3(x,y,z) = 1 -T 2 (x,y) -T 2 (x,z) -T 2 (y,z) + T 3 (x,y,z) 
Pi(x,y,z,t) = 1 -T 2 (x,y) -T 2 (x,z) -T 2 (x,t) -T 2 (y,z) -T 2 (y,t) -T 2 (z,t) 
+T 3 (x, y, z) + T 3 (x, y, t) + T 3 (x, z, t) + T 3 (y, z, t) 
+T 2 (z,y)T 2 (z,t) + T 2 (x,z)T 2 (y,t) + T 2 (x,t)T 2 (y, z) - T 4 (x,y, z,t). 

Substituting these relations into (5.15) and using that L Gn = 0, we obtain (writing the 
terms in the same order) 



(5.16) 



K i^j,ai,aje[0,N] 



[0,7V] 2 



G N (x - y)T 2 (x,y) dxdy 



+ 2 1 Gn(x — y) Gjv (z — t)T 2 (x,z)T 2 (y,t) dxdy dzdt 

'[0,7V] 4 ] 



[0,7V] 4 



Gn{x — y)Gpj(z — t)T4(x, y, z, t) dx dy dz dt 



+ 4 / Gn(x - y)G/v(x - z)T 3 (x , y , z) dx dy dz 

'[0,7V] 3 



-4/ Gn(x - y)G/v(x - z)T 2 (y , z) dx dy dz 

./[0,7V] 4 

+ 2 1 Gn{x - y) 2 dxdy -2 l Gn(x - y) 2 T 2 (x,y) dx dy. 

./[0,7V] 2 ./[0,7V] 2 

Similarly (and as we have seen in the proof of Theorem [TJ we have 

ES^ G N (cii - aj) = - / G N (x -y)T 2 (x,y)dxdy 
77^ J\0.N\ 2 



and the result follows. 



□ 



5.1 The one-dimensional case 

Theorem 4. For the sine-f3 processes with (3 = 1,2,4 as described above, we have 



lim Var(Wiv) = 0. 

7V->oo 
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Proof. Since we already know from Proposition 3.2 that lirn/v-»oo EWjv exists, and in view of 



(2.53), it suffices to show that 



(5.17) 



lim | E 

JV->oo N 2 



Cbi — Qi 



i^oj.Oj-efQ.JV] 



0. 



We apply Lemma 5.1 and now deal with all the terms I\ to /6 in (5.5)-(5.10). First, we 
have 

(5.18) h = 2 G N (x-y) 2 dxdy = 2N 2 (log |2sin7r(x - y)\f dxdy. 

J\0,N] 2 ifO.ll 2 



'[0,N] 

For I 2 , using the explicit expression for Gm, making the change of variables x 1 = x/N,y' 
y/N,t = y — z, and recalling that T 2 is translation-invariant since the process isQwe find 

(5.19) I 2 = -4 G N {x -y)G N {x - z)T 2 (y - z)dxdydz 

J\0,N] 4 

log |2 sin7r(x — y)\ log 



-4A 



[0,1] 2 J[N(y-l),Ny] 



2 sin 7T I x — y + 



N 



T2{t) dx dy dt. 



One may like to think that / f(x - y + jj)T 2 (t) dt ->■ f(x - y) since J T 2 = 1. However, 
log 1 2 sin • I is not regular enough to apply this reasoning and J T 2 may converge only condi- 
tionally. First notice that the cluster functions given in ( 3.16| )-(3.17)-(3.18) satisfy 



(5.20) 



(5.21) 



\T 2 (v)\=0 



Pick two exponents a, b > with a + b < 1. Let us examine the t integral in the right-hand 
side of Q5.19P . Assume first that [-N b ,N b ] C [N(y - l),Ny] and that \x - y\ > N~ a . Note 
that for this to be satisfied it suffices that 

(5.22) (x, y) G S N , a := {(x, y) G R 2 : N~ a < y < 1 - N~ a , \x - y\ > N~ a }. 

By the mean value formula, we may write for some \9\ < N b 



log 



2 sin 7T [ x — y + 



N 



log |2sin7r(x 



7rcos7r(a; — y + 



N ■ 



A^sin7r(x — y + jj) 



0(N 



a-l\ 



since we assumed \x — y\ > N a . Thus 
(5.23) / log |2 sin7r(x — y)\ log 

J\-N b .N b ] 



2 sin 7r I x — y + 



iV 



T 2 {t)dt 



[ log 2 |2sinvr(x - y)|T 2 (t) + (^(Af^ 1 log AT) 

J[-N b ,N b ] 



log |2sin7r(x - y)\(l - 0(N~ b )) + 0(A a_1 log A) = log |2 sin vr(x - y) | + 0(iV" 6 log A) . 



1 hence by abuse of notation we write T2(x,y) —T-i(x — y) as for Theorem [l] 
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where we have used (5.20), then (5.21 ) and |x — y| > N a . We then claim that if |x — y| > N a 
we have 



(5.24) / log |2sin7r(x - y)\ log 

J\t\>N b 



2 sin 7r I x — y + 



T 2 (t) dt = o(l). 



Assuming this, and combining with (5.23) we obtain that 

/ / log |2sin7r(a; — y)\ log 

J(x,y)e\0,l] 2 nS N . a J\N(y-i),Ny] 



2 sin 7T \x — y + 



N 



T 2 (t)dxdy dt 



I, 



(x, y )e[o,i]2nS N>a 



log 2 |2sin7r(x — y)\ dx dy + o(l) 



But it is easy to check, since the integrals converge and (5.20) holds, that the contributions 



of the set where (5.22) does not hold are o(l) as N — > oo. We may thus conclude that 



(5.25) 



-AN 1 



(x,y)e [0,1] 



log 2 |2sin7r(x — y) \ dxdy + o(N 2 ). 



To finish with this I 2 term, it remains to prove (5.24). For /3 = 1,2 this is immediately true 
since T 2 (v) = 0(\v\~ 2 ). For /3 = 4, we notice that the same argument that was used above 
to restrict to \x — y\ > N~ a can be used to restrict to \x — y + jr\ > N~ c (note that the 



initial integral is symmetric in y and z). Inserting the formula for T 2 (3.18), and neglecting 
the 0(l/t 2 ) part of T 2 , we thus have to prove that 



/ 

J t 



d sin2vrt 1 f 2 ^ sint 



\t\>N b ,\x-y+±\>N- c dt 2itt 2ir J t 



log 



2 sin 7r I x — y 



N 



o(l) 



We integrate by parts and find that the boundary terms are negligible, and there remains to 
show that 



/ 



\t\>N\\x-y+U>N- 



sin27rf d 
2vrf di 



log 



2 sin 7r ( x — y + 



1 [ 27Tt sint 



2?r 



dt = of IV 



If the derivative falls on the second factor, we are back to the 0(l/i 2 ) situation which gives 
a negligible term, and for the other term we use 



0_ 

dt 



log 



2 sin 7r I x — y + 



N 



0{N- 



c-l\ 



by explicit computation, which gives that the integral is 0(N C ^ 1 log N) = o(l). This com- 
pletes the treatment of I 2 . 

We turn to I3. Using a similar change of variables, we may write this term 



h = 2N z / / / log|2sin7r(x-y)|log 

'[0,1] 2 J[N(x-l),Nx] J[N(y-l),Ny] 



2 sin 7r I x — y + 



v — u 



N 



x T 2 (u)T 2 (v) dx dy du dv . 
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Very similar manipulations to those above show that log |2sin7r (x — y + can be re- 

placed by log |2sin7r(x — y)\ with a o(N 2 ) correction. This leads us to 

I 3 = 2N 2 / log 2 |2sin7r(x-7/)|(ixdy + o(iV 2 ). 
-'foal 2 



For I±, we have 



/ 

.7 [0,7V] 2 



log 2 



2 sin 



tt(x - y) 



N 



T 2 (x-y) dx dy = o{N 2 ). 



2 sin 



■k(x— y) 



N 



Indeed, we may take away a <5iV-neighborhood of the diagonal, outside of which log 

is bounded by log 2 N and f \T%\ is controlled by log N, using (5.20). Thus the whole integral 
is controlled by iVlog 3 N = o(iV 2 ). 

Adding the above results we find that 



h + h + h + h 



o(N 2 



It remains to show that 1$ and Iq also give o(N 2 ) contributions. The expressions I5 and 
Iq are estimated using explicit formulas for cluster functions of the sine-/3 processes. First 
returning to (3.11 )-(3.12 ) we see that for (3 = 1,2, the entries of and are 0( 1+ ^_^| ) . 

in both these 



Combining with (5.3)-(5.4), it follows that T%{x, y,z) = O [ 1+ \( x 
cases. For (3.13), we have 

K^(x,y) = 



i 



\-+\{x-y)(y-z)(z-x)\ 



o 



(5.26) 

We thus obtain 
T 3 (x,y,z) = O 



i 



l+\x-y\ 



1 



1 + |0- y)(y - z)\ 



0(1) 



+o 



O 

o 



1 



l+\x-y\ _ 
1 

\+\x-y\ 



1 



1 + \{y - z){z - x) 



+o 



1 



1 + \(x - y){z - x)\ 



In ( |5.9| ) we may first (as above) remove a small neighborhood of the diagonals, off of which 
\Gn(x — y)\ and \GN(y — z)\ are bounded by O(logiV). It then remains to estimate 



/ 

J[0,N] 3 



\Ts(x, y, z) \ dx dy dz. 



Replacing T3 by its above estimates, and changing variables to x + y + z and successively 
two out of x — y, y — z, and z — x, we find Jj 7V] 3 l^3( x >y> z )l dxdydz < 0(iVlog 3 N), and 
h = o{N 2 ). 

We finally turn to Iq. The formula for T4 is given by (5.3)-(5.4). Comparing to (3.11)- 



(3.12), we see that for j3 = 1, 2, we have 
T 4 (x,y,z,t) = O 



\ + \{x-y){y- z){z-t){t- z) 
The same reasoning as for I5 gives Iq = o(N 2 ). 
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For P = 4, in the formula for T4 obtained with (3.13), in view of (5.26), there are terms 
which a priori have insufficient decay: they are terms of the form 



d sin7r(xi— X2) d sin 2tt(xs — x±) 
dx\ 2tt(xi — X2) dx-s 2tt(x3 — X4) 



2tt(x2—X3) 



sint 



dt 
t .10 



2tt(x4— xi) 



sint 



dt. 



where (xi, . . . , X4) is a permutation of the variables x, y, z, t. This leads to two different types 
of integrals 



(5.27) 



log 



[0,N] 4 



2 sin 



vr(x - y) 



N 



log 



2 sin 



ir(z - t) 



N 



2ir(x-y) 



sm S 



ds 
s .10 



2n(z-t) 



and 
(5.28) 



log 



[0,N]i 



2 sin 



vr(x - y) 



N 



log- 



sins d sin27r(a; — z) d sin27r(y — t) 

ds ' ^~^w ^^~^w -—dxdydzdt 

s ox 2ir{x — z) ay 2n{y — t) 

tt(z - t) 



2 sin 



N 



2n{x-z) g - ng r2n(y-t) 

ds \ 

s Jo 



sins d sm2n(x — y) d sm27r(z — t) 

ds ■ — -. r — — -. r — dx dy dz dt. 

s ox 2tt{x — y) oz 2-k(z — t) 



For (5.27), we may again restrict the domain to \x — y\ > N a with < a < 1. Then, 



integrating by parts in x gives boundary terms which are negligible, and a new integrand with 



extra decay, involving ^- J 2n ^ x y ^ ^j^- ds or log 
contributions. 



2 sin 



N 



. This leads again to o(N 2 



For (5.28), we may first restrict the integral to \x — z\ > N a and \y — t\ > N a , using 
arguments as above. Then, we may replace Jq 2 "^ ^ ds and Jq™^ ds by ^sgn(x-z) 

and ^sgnijj — t) respectively, making only a o(N 2 ) error. Then note that the integrand in 
x — y is a locally odd function, so we can remove domains \x — y\ < N b , \z — t\ < N b from the 
integration domain. Finally, integration by parts in x gives additional decay, yielding o(N 2 ) 
contribution. We conclude that Iq = o(N 2 ), and the result follows. □ 

Corollary 5.2. For all point processes with finite lim7v-s.oo E(Wjv) and liniTv^oo Var(Wjv) = 
in L 2 (ft) and thus in probability. 



lim EWn as N —?■ 00 

N^oo 



The proof is immediate. Processes satisfying these assumptions and having different 
values for lmiAr^oo EWjv are thus mutually singular, such as /3-sine processes with different 
£€{1,2,4}. 



5.2 The two-dimensional case 



Theorem 5. For the determinantal random point process with kernel (3.32) we have 



lim Var(HV) = 0. 

N—too 
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Proof. The starting point is again Lemma 5.1 We note that in view of (3.32), (3.2) and 
(5.1), all the cluster functions for that process are exponentially decreasing when viewed as 
functions of pairwise distances between arguments. We start with the term (5.6). Replacing 
Gn by n^Ejf, using the translation invariance, and changing variables as before (y — z = it), 
we find 



h + 2Ji 

2N 2 



7T 



[0,1] 2 JNy-[0,N] 2 



E N (N(x-y)) (E N (N(x - y)) - E N (N(x -y + u/N))) T 2 {u)dx dydu. 



Noting that in view of the definition of E (cf. ( |2.35 )-( 2.36| )) we have E^{Nx) = E\(x) and 
from (3.31), E\ behaves like Clog|x| near x = (and similarly near points of the lattice Z 2 ). 
Given rj > there thus exists 5 > such that 



/ / 

J [0,l] 2 n{\x-y-Z 2 \<S} JNy-[0,N] 2 



E\{x — y) {E\{x — y)) — E\(x — y + u/N)) T2(u)dx dy du < rj. 



On the other hand, still in view of the definition ( 2.35 )-( 2.36), E\ is uniformly continuous 
away from Z 2 , hence we may write, as TV — > oo, 

/ Ei(x-y) {Ei{x -y))- Ei(x-y + u/N)) T 2 (u)dx dydu = o(l). 

[0,l] 2 \{\x-y-Z 2 \<6} JNy-[0,N] 2 
Since this is true for any rj, it follows that 

h + 2h = o(N 2 ). 

Similarly, 
h ~ h = 



[0,1] 2 JNx-{0,N] 2 JNy-[0,N] 



E^x-y) [E 1 [x-y + 



u — v 



N 



E\{x — y) I T\(u)T2{v) dx dy du dv . 



The same reasoning shows that this is o(N 2 ). For I4, the change of variables x' = x/N and 
y' = y/N yields 

I 4 = -2N 2 [ E x {x-yfT 2 {N{x-y))dxdy. 

J[0,1] 2 

We may take out a 5 neighborhood of the diagonal and its translates by Z 2 , off of which 
E\{x — y) can be bounded by Clog \x — y\ and T 2 {^jf-) by e~ CN ' 2& ' 2 . The whole term is thus 
o(N 2 ). 

We turn to From Q and ((3732J) we find that \T 3 (x,y,z)\ < e -C{\x-y\ 2 +\y-z\ 2 +\x-z\ 2 ) _ 

As above we have 



h 



[o,i] 3 



E ± (x - y)Ei(x - z)T 2 (Nx, Ny, Nz) dx dydz 



and as above we may take out ^-neighborhoods \x — y\ < 5 or \x — z\ < 5 or \y — z\ < 5 
(and their translates by Z 2 ), outside of which the E\ terms are bounded by log's and T3 by 



-CN 2 S 2 ^ The whole term ig thug ( N 2y 

A very similar reasoning applies to Iq. 
Combining all these, we find the result. 



□ 
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Remark 5.3. In view of the proof above, the same result holds for any process such that the 
cluster functions decay sufficiently fast away from diagonals, for example exponentially. 



6 Miscellaneous computations 

In this section we gather various additional computations of expectations and additional facts. 
6.1 Operations on processes 

In this subsection, we examine the effect on lmi/v-s-oo EWjv of two common operations on 
independent processes: superposition and decimation (see |D V J] ) . 

Proposition 6.1. Let X\,...,Xm be M independent translation invariant point processes 
with density 1 and two-point correlation functions {p2^}k<m> satisfying the assumptions of 
Theorem [JJ 

Assume that limjv^oo EVW(A'j) < oo for i = 1, . . . , M . Let X denote the superposition of 
independent processes Xi, where Xi denotes the image of the process X{ under the dilation by 
factor M of the line. Then, with obvious notation, 

1 M 

lim EW N (X) =logM+ — V lim EW N (Xi). 

N— >oo M N—>oo 

i=l 

(i) -(i) 

Proof. Let T 2 be the second cluster functions corresponding to Xi. Let p 2 now denote the 
second correlation function for the process Xi, which has density 1/M. We have p 2 (x,y) = 



M 



~Pi\fdi m) - ^ ne P rocess % clearly has density 1, and its second correlation function is 



M M 



P^y)= E « J VE^.) = M(M-i)l + lE#&s 

&je[l,M] i=l ' »=1 

We also denote by T<i the corresponding second cluster function. We thus have 



i=l i=l 



M M 2 ^\ 2 \MJJ M 2 ^ 2 \M, 

i=l i=l 

and it easily follows, with a change of variables, that T 2 {v) dv = 1. In addition, 

/OO /'OO 1 
T 2 (v)\0g\27Tv\dv= Y, T 2 ] (^)^g\27Tv\dv 
-00 J — 00 - -J^ 

1 M roo 1 A/ f°° c\ 

= 17 E / T 2 lQ S I 2 ^ MS I ds = l0 § M + Tf E / T 2 W l0 § \ 27TV \ dv - 
M i=l J ~°° k=l J -°° 

The result follows easily using Theorem [T] □ 
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In dimension 2, reproducing the proof, but replacing the dilations by factor M by dilations 
by factor \/M, we obtain instead the result: 



1 1 M 

lim EW N {X) = -logM + — V lim EWjvpQ). 

i=l 

For example, superposing two independent processes with the same second correlation 
function leads to an increase of limjv^oo EWn by log 2 in dimension 1 and \ log 2 in dimension 
2. Superposition can thus be seen as "increasing the disorder". 

We next turn to decimation. 

One can define a 'random decimation' of a process by erasing points at random with 
probability 1/2. The second correlation function then transforms into 

Rz{x,y) = -p 2 (x,y). 

But the space needs to be rescaled by a factor 2 in order to maintain a density one, so the 
correlation function after that is 

P2(x,y) = P2(2x,2y). 

It is clear in view of Theorems[T]and[2]that if a process has finite lim7v-s>oo EWn, its decimation 
will not, since the condition J T 2 = 1 will be destroyed by this operation. 

On the other hand, one can define a 'deterministic decimation' by erasing every even (or 
odd) point of an (ordered) random point configuration followed by rescaling of the space to 
keep the density at 1. While we cannot say anything about this operation in general, one can 
observe what it does in a couple of cases. 

It is known, see e.g. |AGZ| page 66], that the (3 = 2 sine process is the deterministic 
decimation of superposition of two (3 = 1 sine processes, or symbolically 

(sine (3 = 2) = decimation((sine (3 = 1) U (sine (3 = 1)). 

From Proposition |3 . 2| we know that limAr^oo EWV is 1 — 7 for the left-hand side, and Proposi- 
tion 6.1 says that liniAr_ 5 . 00 EWn is 2—7 for (sine (3 = l)U(sine (3 = 1). Thus, the deterministic 
decimation decreased the value of limjv-^oc EWn by 1. 
Similarly, 

(sine j3 = 4) = decimation((sine f3 = 1)), 



and we see that the decimation decreased the value of lim7v-s>oo EWV from 2 — 7 — log 2 to 

3 
2 



3 - 7 - log 2 



6.2 Discrete {3 = 2 sine process 

The (3 = 2 discrete sine process was first obtained in [BOO] as the bulk scaling limit of 
the Plancherel measure for symmetric groups, and it was shown in [BKMMJ to be the a 
universal local scaling limit for a broad family of discrete probabilistic models of random 
matrix type with = 2. The goal of this section is to compute limAr^oo EWn for the suitably 
scaled discrete sine process embedded into the real line. By the construction, this provides 
an interpolation between the case of the perfect lattice, for which Wn = 0, and the case of 
the continuous (3 = 2 sine process treated in the previous sections. 
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Let p E (0, 1). The discrete sine process with density p is a random point process on 
with the correlation functions (k > 1) 



Pk(xi, ... ,x k ) = det 



shi7r p{x% — Xj) 
ir(xi - xj) 



i k 



Proposition 6.2. Embed Z m£o R wo n i— > pre; i/iis turns the discrete sine process of density 
p into a random point process on M with density 1. For i/ie latter process we have 

(6.2) lim EHV = />logp + 2 V f Sm(p?rn) Vlog(27rH- 



7V-s>oo 



Proof. For the calculation, we assume that y is an integer (the same argument however should 
hold without this assumption by examining more carefully error terms). The calculation is 
then can be viewed as a discrete version of that of Theorem [TJ First, by definition (2.53) of 
Wjv, we have 



N 



E P2(»,i)log 



2 sin 



np(i-j) 



N 



+ logiV 



and since P2(i,j) = p 2 



sin wp(i—j) 



we find 



(6.3) EW N = ± 



SU17T p(i — j) 

ir(i-j) 



log 



2 sin 



irp(i-j) 



N 

2 



+ logiV 



P_ 

N 



E fc* 

ifte[l,N/p] 



2 sin 



iV 



Let us first examine the contribution of the last sum. From the knowledge of WV(Z) = 0, we 
know that for K E Z we have 



lim — 



E 



Applying this to K = N/p we find that 



2 sin 



ir(i - j) 



K 



+ log if = 0. 



(6.4) 



El 

N 



E 

i^ie[i,iv/p] 



2 sin 



N/p 



-p\ogN/p + o{\). 



We next turn to the first two terms in (6.3). As in Theorem [TJ we expect only the near 
diagonal terms to contribute, so that 



log 



2 sin 



irp(i - j) 



N 



log 



2irp(i-j) 



N 



log \2Trp(i — j)\ — logiV. 



Inserting this and (6.4) into (6.3) we find 
(6.5) 



1 

N 



E 



smirp(i — j] 
ir(i - j) 



(log |27T/>(« - i)| - log AT) -plog./V/p + logiV + o(l). 
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Figure 2: Numerical evaluation of limAr_> 00 EWat for the discrete sine process 



We first focus on 

. . - log N 
(6.6) 



E 



N ^ \ ir(i- j) 



sin7rp(i — j) 



N/p-i 9 , 
S E (N/p-u) 2s ' m{7Tpu) 



N 



u=l 



ir 2 u 2 



1 AT N /C- 1 1 o 

log iv ^—v 1 — cos Zirpu . 
> s-s r-o(l). 



u=l 



Indeed, we can bound YlfJi u 2sm ^ n / u ^ by h = OQag N/ p) and this multiplied by ^j^- 

is negligible as N — > +oo. The last sum then appears as a Fourier series and can be computed 
explicitly, which leads to 
(6.7) 



log N N {^-\ 1 1 — cos2irpu log A ( tt 2 1 

6~ ~ 12 



E 

u=l 



12tt V - \2irp + 2tt' z ) ) = (p - 1) log A. 



We next turn to 
(6.8) 

1 

N 



E 



• , \ 2 N/p-l 

1 " ■< \ log|27rp(i-i)| = - £ W/.-U) 



V ^ - i) 



(sin TrpuY 



N 



log(2-K pu). 



u=l 



Again, the term containing u can be neglected since it is bounded by Ylu=i h log(27rpu) < 
0( 



o(l). Combining (6.5)-(6.8) and letting A — >■ oo, we finally arrive at (|6.2[). □ 



N ' 



The graph of limjv^^ EWV in (6.2 ) is presented in Fig. 6.2 , it shows a function decreasing 



from 1 — 7 at p = to at p = 1, as expected. 
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